WEAKLY COUPLED BOUND STATES OF PAULI OPERATORS 
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Abstract. We consider the two-dimensional Pauli operator perturbed by a weakly 
coupled, attractive potential. We show that besides the eigenvalues arising from the 
Aharonov-Casher zero modes there are two or one (depending on whether the flux 
of the magnetic field is integer or not) additional eigenvalues for arbitrarily small 
coupling and we calculate their asymptotics in the weak coupling limit. 



1. Introduction and main results 

1.1. Introduction. This paper is concerned with negative eigenvalues of perturba- 
tions of the two-dimensional Pauli operator 

P := (a ■ (— iV + A)) 2 in L 2 {R 2 , C 2 ) . (1.1) 

Here o = (cr 1; cr 2 ) is the pair of the first two Pauli matrices 

A /0 -i 



n> \i o) ' °* = \i 

and A is a real vector potential corresponding to the magnetic field B = curl A. Of 
course, P is non-negative, and it is well-known that the point may be an eigenvalue 
of P. Indeed, the Aharonov-Casher theorem (see, e.g., (6j Thm. 6.5]) asserts that if 
B is, say, bounded with compact support, then the dimension of the kernel of P is 

N := #{m e N : m < |$| - l} , (1.2) 

where 

$ := — / B(x) dx (1.3) 



2?r„ 

is the total flux of B. It is less known that P also has a virtual level at 0. Indeed, 
it was shown by Weidl [T9] that if V is a non-negative, sufficiently regular and not 
identically zero function, then for all sufficiently small a > the perturbed Pauli 
operator P — aV has exactly 

N'-=[ N + l lf$eR \ Z > (14) 
[N + 2 if$eZ 

negative eigenvalues. We express this fact by saying that P has, in addition to its 
eigenvalues, one, resp. two virtual levels at zero. 



© 2009 by the authors. This paper may be reproduced, in its entirety, for non-commercial 
purposes. 



2 



RUPERT L. FRANK, SERGEY MOROZOV, AND SEMJON VUGALTER 



These 'additional' eigenvalues are of physical interest, in particular, since an anoma- 
lous magnetic moment g > 2 corresponds to a perturbation aV = (g — 2)B of the 
Pauli operator. Note that g = 2.0023 for an electron. We refer to [3] and references 
therein for more on this. 

Let Ai(a), . . . , Ajv'(a) denote the N' smallest eigenvalues of the operator P — aV 
in non- decreasing order. Throughout we assume that B is compactly supported and 
radially symmetric. The goal of this paper is to obtain the asymptotic behavior of 
these eigenvalues as a — > 0+. While it follows in a rather straightforward manner that 

\j{a)~ -cja, j = l,...,N, 

as a — > 0+, our main result is that 

Ajv+i(a) ~ -CTv+ia 1 ^ if /i := {|$|} = |$| - N G (0,1) 



and 



\ N+1 (a) ~ - . C f N+l0c , , ln\X N+2 (a)\ — if|$|=JVeN. 

ln(cjv + ia) c N+2 a 



(The result is slightly different in the case $ = and we refer to Theorem 11.31 below 
for the precise statement.) Moreover, we obtain explicit expressions for the coefficients 
Cj > as well as estimates for the remainders in the above asymptotic expressions. 

We note that our result quantifies a paramagnetic effect of the Pauli operator. 
Indeed, while the ground state energy of the Schrodinger operator —A — aV is expo- 
nentially small in a in the weak coupling limit (see [IE]), the addition of an arbitrarily 
small magnetic field with non-zero flux $ leads to a much more negative ground state 
energy of the Pauli operator P — aV which is of the order a 1 /'*' if |$| < 1, a\ lna| _1 
if |$| = 1 and a if |$| > 1. 

The existence of weakly coupled eigenvalues can intuitively be understood as follows. 
Introducing the function 

£(x) := -(27T)- 1 / B(y) \n\x-y\dy, (1.5) 

one easily finds the well-known relation 

/ \a- {-iV + A)i)\ 2 dx = 4 f (e- 2 «|^e^+| 2 + e 2? |^e-^_| 2 ) dx, (1.6) 
Jr 2 Jr 2 

where ip — (4' + ,' l P~), d z = \{d Xl — id X2 ) and dz = \{d xi + id X2 ). This suggests that 
zero modes (i.e., solutions of the equation Pifj = 0) should be of the form (fijj", 0) and 
(0, fijjr ), where 

Q±( x ) : = (27T)- 1 / 2 e^\ Xl ± %x 2 ) k , k E N . (1.7) 

For the sake of definiteness let us assume that $ > 0. If B is radial and compactly 
supported, then £(x) = — $ln \x\ for large \x\ by Newton's theorem and hence all the 
are increasing at infinity and do not belong to L 2 (M 2 ). In contrast, fl^ belongs 



WEAKLY COUPLED BOUND STATES OF PAULI OPERATORS — March 27, 2009 3 

to /^(M 2 ) iffO^A;< ( l ) — 1. Those fl^ give rise to eigenvalues of P — aV which 
disappear linearly in the weak coupling limit. What is more important for us is that 
the functions with $ — 1 ^ k ^ $ are bounded (although they do not belong to 
/^(IR 2 )). It was already observed in [3] that any perturbation by a negative potential 
— V will turn these functions into Z^-eigenfunctions. Our key point is that the weak 
coupling asymptotics of the eigenvalues are determined by the spatial asymptotics of 
the functions . 

There is an enormous literature on weakly coupled eigenvalues and low energy be- 
havior of Schrodinger operators from which we only mention [T6l HI HH EL7J [HI [19] , the 
surveys [51 [15] and the recent papers [H [2j [121 E] • We emphasize that techniques from 
weakly coupled Schrodinger operators have also turned out to be useful in a non-linear 
context [7]. 

All papers on weak coupling asymptotics, which we are aware of, are based on the 
Birman-Schwinger principle and on operator-theoretic arguments. They rely upon 
very detailed knowledge of the unperturbed Green's function which is, of course, ex- 
plicitly known for the Laplacian. It seems very unlikely that such information can be 
obtained in the generality in which we work here. 

Instead we propose a completely different, purely variational approach. Its ad- 
vantage is that it closely follows the above mentioned intuition that weak coupling 
asymptotics are determined by spatial asymptotics of resonance functions, which is 
often obscured in the operator-theoretic approach. It is by no means restricted to the 
problem under consideration and it allows one to recover in a simple way and improve 
upon the results in [16j HI [TUl [HJ [12] concerning the lowest eigenvalue. While the 
ground state energy is of primary physical interest, we should say that we do not see 
in general how to obtain results on excited states with our method. This is where 
in the present situation the radial symmetry of the problem comes in, which reduces 
the problem to ground state problems for half-line operators. We emphasize, however, 
that our techniques are not restricted to one-dimensional problems. 

Acknowledgements. The authors would like to thank T. Weidl for drawing their at- 
tention to this problem and for helpful discussions. R.F. has profitted from discussions 
with T. Ekholm and H. Kovarik. The authors gratefully acknowledge the hospitality 
at Stuttgart University, KTH Stockholm and ESI Vienna, where parts of this work 
were done. This work was partially supported through Deutsche Forschungsgemein- 
schaft's (DFG) grant FR 2664/1-1, US National Science Foundation's grant PHY 06 
52854 (R.F.), DFG grant SI 348/12-2, Engineering and Physical Sciences Research 
Council's grant EP/F029721/1 (S.M.) and DFG grant WE 1964/2 (S.V.). 

1.2. Main results. Let us state the precise conditions on the magnetic field and the 
electric potential. 
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Assumption 1.1. Let B and V be radial, real- valued and measurable functions with 
compact support in M 2 such that 

/ \B\(l + ln_ |x|) dx < oo and / |V| (l + ln_ \x\) dx < oo . (1.8) 
Jr 2 Jr 2 

Here and below t± := max{±t, 0} for a number or a function t. Let A G L2,i oc (K 2 , M 2 ) 
be a vector field with curl A = B (see (13. 2 p for a convenient explicit choice). Under 
Assumption 11.11 the operators P — aV, a G M, are defined through the closure of the 
quadratic forms 

/ + A)?P\ 2 dx-a [ V\ip\ 2 dx , ip G C^°(IR 2 , C 2 ) . 

Jr 2 Jr 2 

Using $, A, A' and fi^ from (jOjl . (TO) . (Ojl and (fUD we define 

u fc := / U|^| 2 dx, if ± $ > and fc = 0, . . . , A' — 1 , (1.9) 

Jr 2 

v$:=[ V\tt±\ 2 dx, if$ = 0, (1.10) 



and mk := j R2 \^k\ 2 dx if ±$ > and k = 0, . . . , A — 1. The main result of this 
paper are the following two theorems concerning the weak coupling asymptotics for 
Pauli operators. 

Theorem 1.2 (Case of noninteger magnetic flux). Let B satisfy Assumption \l.l\ and 
let V be a non-negative, not identically vanishing function satisfying Assumption \l.l[ 
In addition, assume that $ G K\ Z and put \i := { |$| } . Then for all sufficiently small 
a > the operator P — aV has exactly A + 1 negative eigenvalues Ai (a) , . . . , A tv +1 (a), 
and as a —>■ 0+ one has for j — 1, . . . , A 

v j-i ( , , i°( a ) *fj < l $ l " 1 

and 

,.. ■= ( 

r(i-Ai) 

The result for integer flux takes the following form. 

Theorem 1.3 (Case of integer magnetic flux). Let B satisfy Assumption \1.1\ and 
let V be a non-negative, not identically vanishing function satisfying Assumption 
In addition, assume that $ G Z. Then for all sufficiently small a > the operator 
P — aV^ has exactly A + 2 negative eigenvalues Ai(a), . . . , \n+2(®)- ^4s a — > 0+ £/ie 
eigenvalues Ai(a), . . . , A at (a) satisfy (11. lip and 

Aj-(a) = — ^-a(l + 0(a|lna|)) i/j = |$| — 1 G {1,..., A} . 



A ( a ) = L^ a i i + J ^ ' (1.11) 

J mj_i V l0K) i/j > |$| -1 ' 



A w+1 (a) = -c^a^l + 0(a min{1 ^- 1} )) , 
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Moreover, i/$6Z \ {0} ; then 



a ( ^/ln|lna|\\ 
Ajv+2(a) = -exp ( (l + 0(a) 



2 

v N+1 o 

and if $ = ; t/ien 

2 



Ai(cn) = - exp C- -,— (1 + 0(a)) 

V t) a 

/ 2 

A 2 (a) = - exp - + 0( a ) 

V t> a 



We emphasize that the meaning of the asymptotics in the case $eZ\{0}is that 
In | A («) | = — ^ v N+i a l (l + 0(a)), and similarly in the case $ = 0. 

Remark 1.4. The assumption that V is non- negative can be somewhat relaxed. Indeed, 
our proof shows that for any k G {0, . . . , N' — 1} for which v k > (with the obvious 
modification for $ = 0) the operators P — aV have a negative eigenvalue with the 
asymptotics given in the theorem. This is in agreement with a result of Weidl [19] 
who has shown that for V ^ one has 

lim N(P - aV) = # {k 6 {0, . . . , JV' - 1} : v k ^ 0} 

(with the obvious modification for $ = 0). Here N(P — aV) denotes the number of 
negative eigenvalues, counting multiplicities, of P — aV. We have not been able to 
compute the precise asymptotics of the eigenvalues corresponding to v k = 0. 

Let us comment on our assumptions. We assume that V has compact support 
mainly for the sake of simplicity in order to avoid additional technicalities. With 
additional work one can probably also replace the support assumption on B by a 
suitable short-range decay assumption. The assumption that B and V are radial, 
however, is crucial for us, at least if |$| ^ 1, since it allows us to avoid orthogonality 
conditions in the case of several eigenvalues and instead to work with a finite number 
of ground state energies. While physically reasonable we would find it mathematically 
desirable to remove this assumption. 

1.3. Weak coupling asymptotics for half-line operators. Because of the as- 
sumed radial symmetry of B and V Theorems 11.21 and 11.31 reduce to statements about 
families of half-line operators. In this subsection we shall formulate a weak coupling 
result for a more general class of one-dimensional operators. Throughout we work 
under 
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Assumption 1.5. (a) The real-valued function W G Lx,i O c(0, oo) Q satisfies 

W(r) = fj 2 r~ 2 , r^R, (1.13) 

for some R > and fi ^ 0. 

(b) For some ^ a < 1 and some M ^ one has for all ip G C^°(0, oo) 

/*oo /*oo roc 

/ W_\ip\ 2 r dr ^ a / |^'| 2 r dr + M \ip\ 2 r dr . (1.14) 
Jo Jo Jo 

(c) For all ip G C^O, oo) one has 

(\ifj'\ 2 + W\ip\ 2 ) r dr ^0. (1.15) 



The closure t of the quadratic form on the left hand side of f II . 1 5 j) in the Hilbert 

space L 2 (M + ,rdr) has domain 

domt = {tp e L 2 (M+,rdr) : if;', W]/ 2 ip G L 2 (M+,rdr)} . 

It generates a non-negative self-adjoint operator T which acts on functions ip G domT 
according to 

Til) = -r^ 1 (r^ / )' + Wip. 

We denote by V the space of all functions V G Li i oc (0, oo) with compact support which 
are form compact with respect to T, that is, the operator (T + 1)~ 1 ^ 2 V(T + l) -1 / 2 is 
compact. For V G V the operator T — V is defined as usual via its quadratic form 
with domain domt. By Weyl's theorem its negative spectrum (if non-empty) consists 
of discrete eigenvalues of finite multiplicites. We are interested in the behavior of the 
lowest eigenvalue of the operators T — aV in the weak coupling limit a — > 0. 

Along with the form domain dom t of T we shall use the local form domain of T 
defined by 

domi oc t := |^ G L 2 ,ioc([0, oo),rdr) : ip', W+ 2 ip G L 2 ,i O c([0, oo), r dr) j . 

We shall say that ip is a weak solution of the equation Tip — if ip G domi oc t and 

t[tp, ip\ = for any G domt with compact support in [0, oo) . 

Note that this notion incorporates a boundary condition at zero but no decay condition 
at infinity. 

Here is the key for proving Theorems 11.21 and 11.31 



■^i,ioc(0, oo) denotes the space of functions which are integrable on all compact subsets of (0,oo). 
In view of later applications we insist on the fact that W may have a non-integrable singularity at 
the origin. 
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Theorem 1.6 (Weak coupling asymptotics for half-line operators). Suppose that As- 
sumption 1 1.5\ holds and that 

there exists a positive weak solution tp of the equation Tip = such that 

V> (r)=r~'S r^R, (1.16) 
with jj, and R from (11.131) . 



Let V G V such that 



v :- 



Then for all sufficiently small a > 0, T 
and as a — > 0+ one has 



V(r)^(r)rdr > 0. (1.17) 
aV has a unique negative eigenvalue X a , 



( ( 2 

- exp (l + O 

V , av 



a 



— c M (av) » (l + 0(a 



min{l ,1-1} 



)) 



lav 
I In a\ 



i + o( 



In 


In a 




In a 





aw 



j™i)l(r)rdr 



0(a min{1 ^- 1} ] 
0(a\ lnal) , 



^2 
\i = 2 



// = 0, 
p 6(0,1) 

A* = 1 > 
// > 1, 



where c^ is given by (11.121) . 



We remark that for // > 1 the function ip Q is square-integrable, and hence is 
an eigenvalue of T. In this case, the leading order asymptotics of A Q follow from 
the abstract arguments of [T7]. Since in our approach there is hardly any difference 
between the cases /i > 1 and ^ \i ^ 1, we include an independent proof which, 
moreover, yields a remainder estimate. We note in passing that by concavity the 
O-term is positive for /x > 1. 

In general, the coefficient v in (II. 17ft depends implicitly on the background potential 
W through the function ipQ. When applying Theorem II .61 in the proof of Theorems 11.21 
and 11.31 however, W will have a specific form which allows to determine ipo explicitly. 
Hence for given V the coefficient v can be computed explicitly. 



The main point in Theorem [L6] is that it connects the existence of a positive solution 
of the equation Tipo = with a certain behavior at infinity to the spectral information 
about the existence of negative eigenvalues. This connection is quantitative in the 
sense that the decay of ipo determines the size of the eigenvalue in the weak coupling 
limit. An initial step in our proof of Theorem II .61 will be a qualitative version of this 
correspondence, which we single out as Proposition 11.71 below. We emphasize that 
this qualitative statement is a well-studied feature of second order elliptic operators 
[15] . though our assumptions on the potentials seem to be much weaker than those 
typically imposed in the literature. 
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Proposition 1.7 (Characterization of criticality) . Suppose that Assumption \1.5\ holds. 
Then Assumption (11.161) is equivalent to each of the following: 

(i) For any non-negative ^ V G V the spectrum of T — V in (— oo, 0) is 
nonempty. 

(ii) For some V G V and any a > t/ie spectrum of T — aV in (— oo,0) is 
nonempty. 

The structure of the proofs of Theorem 11.61 and Proposition 11.71 is as follows. In 
Section [27T1 we assume the existence of a function ip as in (I1.16p . Starting from this 
function we construct a family of trial functions with negative energy. This implies 
part (i) in Proposition 11.71 and gives us the upper bound claimed in Theorem 11.61 In 
Section 12.21 we assume the existence of a potential V as in part (ii) of Proposition 
11.71 Starting from the family of eigenfunctions ip a of the corresponding operators we 
construct a function ipo as in (11. 16ft . Moreover, controlling the convergence of ip a to 
ipo will allow us to prove the lower bound claimed in Theorem 11.61 Uniqueness of the 
eigenvalue is the content of Lemma 12.51 

2. Proof of Theorem 11.61 and Proposition 11.71 

2.1. The upper bound. In this subsection we prove that (I1.16P implies (i) in Propo- 
sition 11.71 and we derive the upper bound on \ a stated in Theorem 11.61 Both follow 
immediately from 

Proposition 2.1. Assume that there exists a function tp Q as in (I1.16P and let V G V 
satisfy (11.171) . Then X a := inf spec(T — aV) < for any a > and, as a — > 0, 

- exp ( (1 + const a) V a = , 

\ i av / 

-c M (av)* (1 — const a mm ^ 1 '^~ 1 ^) , // G (0, 1) , 

2av . . 

const — p: j— I , /i = 1 



In a 



av 



[i > 1 



I So i>l{r)rdr' 
where c^ is given by (11.121) . 

Proof. First assume that // > 1. Then ip is square-integrable and t[ipo] = 0. Hence 
by the variational principle 



t[ip ] ~ a So°° V{r)4>l{r)r dr _ , , , 



K ^ f oo , 2( v , = - / ^ ( r ) rdr aV > 

Jo Vw) r dr V ^0 7 

as claimed. In the remainder of this proof we shall assume that ^ /x ^ 1. Let R > 
such that W(r) = /i 2 r~ 2 and V(r) = for r ^ i?/2 and define for any k > 



^o(r) , r < R 



<P*{r)-=< K^kt) (2.1) 
" r > R. 
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Here is the modified Bessel function of order //, see [8]. The function ip K belongs 
to the form domain of T — aV and we claim that for small k > 



t[<p K ] 



-(Inn)- 1 + 0(|ln/?|- 2 ) 



and 



Wk\ 



0(k 2 



0(k- 2 \\tik\- 2 ) , 
2 i- 2 ,/(l " *0 



Inn + 0(1) 



K 



2(1—/*) _|_ Q^2min{-l+2 l u,0}^ 



/i = 0, 

^e(o,i) 
/" = l, 

fie (0,1) 

/U = 1 . 



(2.2) 



(2.3) 



To prove the first relation, we multiply the equation Tijj Q = by rip Q and integrate 
by parts. Using the definition of a weak solution and that ipo e domi oc t one sees that 
there appear no boundary terms at zero and one obtains 



Hence 
with 



/ (l^o| 2 + ^o| 2 ) rdr = -(lFT 2 ». 
Jo 

t[<p K ] = -fiR- 2 ^ + RT 2 »K- 2 {kR) A{k) 



(2.4) 
(2.5) 



From [8, 5.54.2, 8.486.13] we get 



2 (K,r) )r dr . 



A(k) = / (K 2 rK 2 {Kr) - 2 K (iK fl+1 {Kr)K fl {Kr) + 2^ArJ(/cr)) d; 

/ K,rK 2 +1 (Kr) d(nr) + 2(i / K ^(kt) d K ^(kt) 
Jr Jr 

■(K,(kR)K, +2 (kR) - KI +1 {kR)) - (iK 2 (kR) , 



iR 

K 2 R 2 



and ( 12. 2D follows from the asymptotics of the modified Bessel functions for small 
arguments % 8.485, 8.445, 8.446]. To prove ( T273T) we write 

■R 



[ \ifj \ 2 rdr + R- 2 ^K- 2 {KR)J{K) 
Jo 



where 



J(«) 



( KI{kt)v dr = —(K^kR^^kR) - K^kR)) . 



In the last identity we used again [H 5.54.2], and asymptotics ( 12. 3 j) follow as before 
from the asymptotics of the modified Bessel functions. 



K 2 
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Relation ( 12.2ft shows that for any a we can choose a sufficiently small k to make the 
quotient (t[ip K ] — av)/\\<p K \\ 2 negative. To obtain an upper bound on A a we minimize 
this quotient to leading order in a. Namely, we choose 

'Cexp(-2/(W)) , fi = 0, 

' /iG(0 ' 1) ' 

2av 

— ; , = 1, 

k — In a 

where C is a sufficiently small constant in case /i = 0. The claimed upper bounds 
easily follow using (E2D and (EJ). □ 

Remark 2.2. The definition of </2 K also makes sense for /x > 1. For later reference we 
note that in this case ( 12. 2ft and ( 12.31) take the form 

L J \0( K 4 |ln K |), /i = 2, V ^ 

and 

m K 2 = ^0 + < 2.7 
This is proved similarly as in the case ^ /i ^ 1. 

2.2. The lower bound. Our goal in this section is to prove that (ii) in Proposition 
11.71 implies the existence of a function ipo as in ( 11.161) and to prove the lower bound 
stated in Theorem 11.61 

Throughout we assume that T is non-negative and we fix V G V such that the 
negative spectrum of T — aV is non-empty for any a > 0. By Weyl's theorem, 
A a := inf spec(T — aV) < is an eigenvalue and we normalize the corresponding 
eigenfunction ip a by 

ip a (R) = R-v . (2.8) 
Here R > is chosen such that W(r) = fi 2 r~ 2 and V(r) = for r ^ R/2. 

2.2.1. Existence of a virtual ground state. We prove that ip a , normalized by (12. 8p . 
have a limit as a — > +0 and that this limit is a weak solution of the equation Tip = 0. 
More precisely, we have 

Lemma 2.3. The functions ip a converge pointwise and in I/ 2i i O c([0, oo), r dr) to a 
function ip satisfying the properties stated in (11.161) . Moreover, 

r R 

0</ (\(ip a -ipo)'\ 2 + W\ip a -ip \ 2 ) rdr ^ const a 2 (2.9) 
Jo 
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for all sufficiently small a > 0, and 



R r R 

2 „ a „ _ / U. 12 



\vjj a \ 2 rdr= / |V» rrdr + O(a), (2.10) 

R pR 

V\tp a \ 2 r dr = / l/|Vo| 2 rdr + 0(a). (2.11) 
Jo 

For the proof of this lemma we need a rough a priori lower bound on the lowest 
eigenvalue, which we single out as 

Lemma 2.4. For all sufficiently small a > one has 

\X a \ ^ const a . (2-12) 

Proof. First note that X a — > as a — > 0. Indeed, this follows from the fact that the 
operator {T + e)~ l l 2 V (T + e)~ l l 2 is bounded for any e > 0, which in turn follows from 
the form boundedness of V. Now by the variational principle, X a is the infimum over 
linear functions of a, and hence X a is a concave function of a. Hence X a ^ aX% for 
< a < 1, proving (12121) . □ 

Proof of Lemma \2. 31 First step. Definition ofipo. Let Tr be the self-adjoint operator 
in ^2(0, .R, r dr) associated with the closure of the quadratic form 

*i#]:= r(l^'| 2 + ^|^| 2 ) r dr, 1>eC?(0,R). 
Jo 

Using (11.141) one easily shows that T R has discrete spectrum. By (11.151) T R is non- 
negative. We claim that Tr is actually positive definite. Indeed, otherwise T R would 
have a zero eigenvalue with corresponding eigenfunction tp . The extension of ipQ by 
zero belongs to the form domain dom t and according to (11.151) minimizes the quadratic 
form t. Hence it belongs even to the operator domain domT. But this would imply 
that i)'o(R) = and therefore by the unique solvability of the Cauchy problem ip = 0, 
a contradiction. 

Fix a function F G domT with F(R) = R~^ and denote / := TF. We define on 
(0,R) 

:=F-T~ l f. 

(Strictly speaking, T^ 1 is applied to the restriction of / to the interval (0, R).) Note 
that ipo satisfies 

-r-\r^)' + H/^ = in (0, R) , fo(R) = Br" , (2.13) 

and belongs to domi oc )f:ij (which is defined similarly to domi oc t). Moreover we claim 
that, 

rR 

/ (W \ 2 + W \M 2 ) r & r is finite. (2.14) 
Jo 
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Indeed, 

r R 

/ {Wo\ 2 + W\M 2 ) rdr 
Jo 

= t R [ip — F] — I (\F'\ 2 + W\F\ 2 ) r dr + 2 Re / (F 7 ^ + WF^ ) r dr 
Jo Jo 

rR_ pR_ 

= t R [ip - F] - / fFr dr - RF'(R)F(R) + 2 Re / f^ r dr 
Jo Jo 

+ 2RReF'(R)tfjo(R) , 

where all the terms on the right hand side are finite. (Indeed, one easily checks 
that F'(R) is well defined for F e domT.) Note that when integrating by parts no 
boundary terms at zero appear since F and ip belong to domi oc t R . 

Second step. Proof of (\2.9\i . Let ^ x ^ 1 be a smooth function on [0, R] such 
that x = 1 on [0, R/2] and x = on [3R/4, R] . We first claim that 

^ Q -^o = <Pa := T R 1/2 {a{T R 1/2 VT R 1/2 )T 1 R /2 X ^ a + X a T^ a ). (2.15) 

Note that this identity shows that the definition of ipo is independent of the choice of 
F. On the other hand, the above formula is also independent of the choice of x- 

Since T R l ^ 2 VT R 1 ^ 2 is a bounded operator and x^a belongs to the form domain of 
T R , the function ip a is well-defined and belongs to the form domain of T R . For any 
ip G domt R one has (with <p denoting its extension by zero) 

f R 

t R [<p,(p a ]= / (aV + X a )<f4) a r dr = t[<p, ip a ] = t R [(p,il) a - F) + (<£,/) 



= t R [ip, ip a - F}+ t R [ip, F — ipo] = t R [tp, ip a - ipo] . 

Since T R is positive definite this implies (p a = ip a — ip , proving (12.151) . 

We denote by || • \\ R the norm in L 2 (0, R,rdr). It follows from (I2.15p . (I2.12j) . the 
form-boundedness of V and the positive definiteness of T R that 

t R [ip a - M 1/2 = *rM 1/2 < const a(\\i/> a \\ R + t R [ X ^a} 1/2 ) ■ (2.16) 

To estimate the first term on the right hand side of (12.161) we use that T R is positive 
definite and hence 

UJr < UoWr + hah < const (1 + t R [ip a f 2 ) . (2.17) 

To estimate the second term we use that according to (12.141) (recall that W(r) = fi 2 r~ 2 
if x (r) < 1) 

«C const (J (\^' \ 2 + W\^ \ 2 ) rdr + \\1> \\ 2 R + t R [<p a ] + ||^|||) 

«C const (l + t R [<p a ]) . (2.18) 
Combining (12161) . fl2TT7]) and fl2TT8|) we obtain fT2T9|) . 
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Note that again by positive definiteness, (12.91) implies that 

-R 

— ^o| 2 r dr ^ const a 2 , (2-19) 



o 

which implies (I2.10p . Relation (12.111) follows from (12.91) by writing 

-R pR pR pR 

V\tjj a \ 2 r dr — / V|^>o| rdr= / V\ip a \ 2 rdr + 2 / VipaX^ordr 
o Jo Jo Jo 

and using that T R l ^ 2 VT R 1 ^ 2 is bounded and that x^o belongs to the form domain of 

Tr. 

Third step. Pointwise convergence. By explicit solution, 



Recall that ipo satisfies (I2.13P and hence belongs to the two-dimensional space of solu- 
tions of —ipQ+ii 2 r~ 2 ipo = on [R/2, R\. From convergence (12 . 1 9[) and the asymptotics 
of the Bessel functions we conclude that ipo( r ) — r il 011 [-R/2, -R] and we can extend 
■00 to (R, oo). 

Recalling fl 2 . 1 3 [) and (12.121) we see that the coefficients in the differential equation 
satisfied by ip a converge in L l loc as a — > to those of the equation satisfied by 
ipQ. Moreover, for any fixed r ^ R/2 one has ip a {r) — > ipo{ r ) and ip' a (r) —> ip' (r). 
By standard ODE results (see, e.g., [20, Thm. 2.1]) this implies that ip a converge 
pointwise to ipo on (0,r). 

Finally, ip a are the eigenf unctions corresponding to the lowest eigenvalue and there- 
fore non-negative. Hence their pointwise limit ipo is so as well. By standard ODE 
results (an elementary Harnack's inequality), ipo is positive on (0, oo). □ 

Lemma 2.5. For all sufficiently small a > the operator T — aV has only one 
eigenvalue. 

Proof. Imposing a Dirichlet boundary condition is a rank one perturbation of the 
resolvent and can create at most one negative eigenvalue. Hence it suffices to prove 
that the operator T — aV with an additional Dirichlet boundary condition at R is 
non-negative. This is obvious for the part on (R, oo). The part on (0,R) coincides 
with the operator Tr — aV from the previous proof. Since Tr is positive definite and 
T R X I 2 VT R X I 2 is bounded, we obtain the claim. □ 

2.2.2. The lower bound. In the proof of the upper bound we have used identity (12. 4p for 
the virtual ground state ipo- For the proof of the lower bound we need a corresponding 
inequality for all ip. This is the content of 

Lemma 2.6. For any if> G domt 

'\ 2 + W\ip\ 2 )r dr ^ -fi\tP(R)\ 2 . 
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Here we use R as denned at the beginning of this section, but any R with W(r) = 
(yu/r) 2 for r ^ R would do. Note also that the value if)(R) is well-defined by the 
embedding theorem. 

Proof. If ip(R) = 0, the assertion follows since T ^ 0. Hence by homogeneity, we may 
assume that ip(R) = R~^. For e > we define ip £ (r) := ip(r) for ^ r ^ R and 

^(r):=r-"e- e(r_Jl) I r ^ R. 

Since T is non-negative and ip e G domi, we get 

pR poo 

/ (|^'| 2 + ^|^| 2 )r dr ^ - / (\ij' £ \ 2 + fi 2 r~ 2 \^ £ \ 2 )r dr =: -I £ . 

J0 J R 

Elementary calculations show that I e — » fiR~ 2fl as e — > 0, which proves the assertion. 

□ 

Now we are ready to complete the proof of Theorem 11.61 
Proposition 2.7. Lei 7 G V satisfy (11.171) . Then X a := inf spec(T — aV) satisfies 
- exp ( (1 — const a) ) , ji = , 



i 



-c M (at;)' 1 (1 + const a™' 1 ''' ^) , /x G (0, 1) 

'lav . . 

1 + const — 7- j— ) , /i = 1 



A a ^ < 



In a 



Of / I a mm{l^-l} ; , 

1 + const < , /i > 1 



Jo ^o( r ) rdr \ |a|lna|, y, = 2 

where c M zs jiuen fry (11.121) . 

Proof of Proposition \2. 7[ As before, let tp a be the eigenfunctions corresponding to \ a 
normalized by (12.81) and let (p K be the functions defined in (12.11) . Note that 

Mr) = <P* a (r), r^R/2, (2.20) 

where X a = — n 2 a . In order to find a lower bound on \ a we write 

POO 

KUaW 2 = t[i; a ] - a V\^ a \ 2 rdr. (2.21) 



Using Lemma 12.61 ( 12.4ft and ( 12. lip the right hand side can be estimated from below 
according to 



t[ip a ]-a V\ip a \ 2 rdr > t[(p K J - a / V^jjr dr - const a 2 , (2.22) 
Jo Jo 

and in order to estimate the left hand side we use (I2.10I) and obtain 

Ua\\ 2 > H^JI 2 - const a. (2.23) 
Plugging (12T22D and (12T231 into (l2T2T|) and using the a-priori bound (12TT21 yields 

fcj - A a ||y? K J| 2 ^ aw (1 + const a) . (2.24) 



WEAKLY COUPLED BOUND STATES OF PAULI OPERATORS — March 27, 2009 



15 



The assertion will now be an easy consequence of the behavior of t[</?«J and H^kJI 2 
which was established in the previous section. 

Indeed, assume first that < \i < 1. Using (12. 2j) . (12. 3ft and the fact that a ^ 
const | Adl^, which follows from Proposition 12.11 we deduce from ( 12.241) that 

2 i-2 M r (l^—^M) - const |A Q | min{1 -^ } ) < av . 

Together with the a-priori fact that A a — > 0, established in Lemma 12.41 one easily 
obtains the assertion in the case < // < 1. 

In the cases p = and /j ^ 1 we proceed similarly and we only sketch the necessary 
changes. In order to remove the a 2 -term from the right hand side of (12.241) we use 
the rough bounds a ^ const |ln|A a || 1 if p = 0, a ^ const |A a || In |A a || if p = 1 
and a ^ const |A Q | if p > 1, which are deduced from Proposition 12. 1[ Moreover, if 
/i = 0we estimate ||y? K J| 2 ^ in (12.241) and if /i ^ 1 we estimate £[<£>«J ^ 0. We use 
asymptotics (12.21) for p = 0, (12. 3D for p = 1 and (12. 7h for p > 1. Finally, in case p = 1 
to pass from the bound |A a ln|A a || ^ 2av(l + const (lnal" 1 ) to the claimed lower 
bound we use the following Lemma [2.81 □ 

In the previous proof we used 

Lemma 2.8. Let f(t) = -tint for < t < e _1 and f- 1 the inverse function. Then 
f^ 1 is increasing and 

1 / \ —s ( „ /In I In s\ \ \ 

r W = i^( 1+ °(w)) flSS ^ 0+ ' 

The proof of this lemma is elementary and will be omitted. 

2.3. Absence of negative eigenvalues. In this subsection we briefly comment on 
the case where Assumption 1 1 . 51 is satisfied but the equivalent conditions in Proposition 
II .71 fail. We will need this in the proof of Theorems 11.21 and II .31 Note that if the weak 
solution of Tip = is different from r _At on [R, oo), then it has to increase like lnr if 
p = or like r M if p > 0. 

Proposition 2.9. Assume that there exists a positive weak solution ip of the equation 
TipQ = such that the limit lim T ,^ oo (lnr)~ 1 ^ (r) if p = and lim^oo r _M ^ ( r ) iff 1 > 
exists and is non-zero. Then for all ip G Co°(0, oo) and for all non-negative measurable 
V one has 

\ , r r r 00 \ — i r°° 

#]^- sup / V^pdp ^ 2 p~ l dp) / V\i>\ 2 rdr. 

4\0<KooJo Jr J JO 

Note that the integral f °° ^ 2 p~ l d p is finite for any r > because of the assumed 
growth of ipQ. The above supremum, however, may be finite or infinite depending on 
V . In the latter case one actually can show, arguing as below, that there is no c > 
such that t[ip] ^ c J °° V^|-^| 2 r dr for all ip. 
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Proof. Writing ip = if>Qip and using the equation for ipo we find that 

poo poo poo 

t[ip] = / \ip'\ 2 ip 2 r dr , / V\ip\ 2 rdr = / V\(p\ 2 ip 2 r dr . 
Jo Jo Jo 

The assertion now follows from a classical theorem by Muckenhoupt [HJ Thm. 1.3.1/3]. 

□ 

3. Proof of Theorems 11.21 and 11.31 

3.1. Angular momentum decomposition. In this subsection we exploit the fact 
that B and V are radially symmetric, so that P — aV can be decomposed according to 
the eigenvalues of the angular momentum operator. We introduce polar coordinates 
(r, ip) in M 2 and write 

This establishes a unitary equivalence between ip e L2(M 2 ,C 2 ) and sequences 

(VO e ^©L 2 (M + ,rdr,C 2 ). 



in 



By gauge invariance we may choose the magnetic vector potential A as 

l r 

A(r, (p) := 6(r)(— sin ip, cosy?) , 6(r) := - / B{p)pdp. (3.2) 

r Jo 

Note that A e L 2 ,ioc(K 2 , K 2 ) in view of ( 11.8H . Plugging decomposition (13 .ip into the 
quadratic form of the Pauli operator we find that 

f \ a . (_, V + A)ip\ 2 dx = J2 KM + » 

where 

/■oo 

*£[/]:=/ (l/'r + W^I/Hrdr 
Jo 

and 

W±( r ) := (6(r) + mr" 1 ) 2 ± £(r) . 

We denote by the self-adjoint operator in rdr) corresponding to the form 

We note that Cq°(0, oo) is a form core for this operator. Indeed, by the arguments 
of [18] C£°(R 2 \ {0}) is a form core of (D + A) 2 and since B is relatively form compact 
(see the following subsection), the same is true for P. 
We introduce the functions 

r m e -£(r) if m ^ , 

r m e-« (r) J r e 2€(p) p- 2 ™- 1 dp if m < -1 , 



and 



' r -m gC(r) |J- e -2^(p) p 2m-l dp [f m ^ 1 ? 
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Here £ is the function from ( 11. 51) . which is radial and which by Newton's theorem [T3(, 
Thm. 9.2] can be rewritten as 



/V POO 
B(p) p d p lnr — / B(p)p\npdp, r>0. 



(3.3) 



The functions are important since they solve 

-r~>(u;±)9' + iy±u;± = (3.4) 

and belong locally to the form domain of T*. Moreover, using that £(r) = — $ In r for 
r ^ R and that |£(r)| is bounded for ^ r ^ i?, and assuming $ ^ for the sake of 
defmiteness one easily finds the two-sided bounds (with constants independent of m 
and r) 



m > 
m < 



m 1 r m (r + i?) 
r H( r + i2)-* > 

r m (r + i?)* , 

-i r M( r + i?)*~ 2 H 



(3.5) 



|m 
|m 
|m 



-i 



H 



,r + -R 
-i r H(r + i?)-* 



(l + ln + (r/i?)) 



m ^ , 

-$ < m < 

m = — $ < 

m < -$ . 



(3.6) 



3.2. Proof of Theorems 11.21 and 11.31 Because of ( 11.81) B± and V are form compact 
with respect to —A and hence by the diamagnetic inequality also with respect to 
(D + A) 2 . (A proof of this fact may be based on Proposition 12.91 and the Sobolev 
embedding theorem.) Hence (W^)_ is form compact with respect to — r _1 9 r r9 r and 
V is form compact with respect to T^, so Assumption 11.51 is satisfied. 

First assume that $ > 0. We claim that for all sufficiently small a > the operators 
T~—aV with — $ ^ m ^ have a unique negative eigenvalue \(T m —aV). Indeed, this 
follows from Proposition 11.71 since for these values of m the functions tu" are positive, 
decay like r _ * _m , belong locally to the form domain and satisfy ( 13.41) . Moreover, 
putting 

POO 

■= / V(r)u±(r) 2 r dr 
Jo 

Theorem 1 1 . 6 1 yields the following asymptotic behavior as«-> 0+. If — $ + 1 < m ^ 0, 
then 

/ 



A(T" - aV) 



f °° uj 2 r dr 
JO m 



1 



0(a) , 
0(a\ In a |) 

OK) , 



-$ + 2 < m ^ 
m = -$ + 2 ^ 
— $ + 1 < m < — $ 
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If m = — $ 



1 < 0, then 



aF) 



I In a I 



l + O 



In 


In a 




In a 





18 



RUPERT L. FRANK, SERGEY MOROZOV, AND SEMJON VUGALTER 



If — $ < m < — $ + 1, then with c M from (11.121) 

A(T- - aV) = -c, (at;-)* (1 + 0(a mta ^i- 1 >)) 

If m = — $, then 

2 



A(T" - «V) = -exp ( - — -(l + 0(a)) 



These asymptotics coincide with those claimed in Theorems 11.21 and 11.31 since for 
k G No one has 

and v k = v^ k if ±$ > 0. 

In order to complete the proof of Theorems 11.21 and 11.31 in the case $ > we need to 
show that there exists an a c > such that for all < a ^ a c the operators T m — aV 
with m < — $ and m > 0, as well as the operators T+ — aV, m £ Z, are non-negative. 
Note that Proposition 11.71 shows that these operators have no negative eigenvalues for 
a > small (since is unbounded), but it does not imply Theorems 11.21 and 11.31 
since it gives no uniformity in m. 

Instead, we will deduce the assertion from Proposition 12.91 by showing that 



>3C 



sup sup / V(u+) 2 pdp / (w+) 2 p 1 dp<oo 

mGZO<r<ooJo Jr 

and similarly with u;+ replaced by uo m and the supremum restricted to m < — $ and 
m > 0. According to the two-sided estimates (13.51) and (13.61) on this is equivalent 
to showing that 



sup sup k r I Vp dp < oo and sup | In r | / Vpdp<oo. 

k& 0<r^R Jo 0<r^R Jo 

These estimates are easily deduced from fll.8p . This completes the proof in the case 
$ > 0. 

The proof in the case $ = is similar, but now both cUq and are bounded 
positive solutions which locally belong to the form domain. The non-negativity for 
m^O follows again by Proposition 12.91 

The result for $ < follows from that for $ > since complex conjugation is 
an anti-unitary operator which switches the sign of B. This completes the proof of 
Theorems 11.21 and 11.31 
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